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ABSTRACT 



With HST FGS 3 we have determined a parallax for the white dwarf - M dwarf 
interacting binary, Feige 24. The white dwarf (DA) component has an effective temperature, 
T e ff ~ 56, 000K. A weighted average with past parallax determinations (ir a b s = 14.6 ± 0.4 
milliseconds of arc) narrows the range of possible radius values, compared to past estimates. 
We obtain Rda = 0.0185 ± 0.0008i?Q with uncertainty in the temperature and bolometric 
correction the dominant contributors to the error. FGS photometry provides a light curve 
entirely consistent with reflection effects. A recently refined model Mass-Luminosity Relation 
(Baraffe et al. 1998) for low mass stars provides a mass estimate for the M dwarf companion, 
M-dM = 0.37 ± 0.20Mq, where the mass range is due to metallicity and age uncertainties. 
Radial velocities from Vennes and Thorstenscn (1994) provide a mass ratio from which we 
obtain M DA = 0.49±°;<^ M.q. Independently, our radius and recent log g determinations yield 
0.44 < Mda < 0.47A4q. In each case the minimum DA mass is that derived by Vennes & 
Thorstensen from their radial velocities and Keplerian circular orbits with i < 90 deg. Locating 
Feige 24 on an M — R plane suggests a carbon core. 

Our radius and these mass estimates yield a -f gra v inconsistent with that derived by Vennes 
& Thorstensen. We speculate on the nature of a third component whose existence would resolve 
the discrepancy. 



Subject headings: astrometry — stars: individual (Feige 24) — stars: distances — stars: binary 
— stars: white dwarfs — stars: late-type 
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Introduction 



Feige 24 ( = PG 0232+035 = HIP 12031) is DA white dwarf, red dwarf (M1-2V) (Liebert fc Margon 



1977) binary (P = 4.23 days, Vermes & Thorstensen 1994=VT94) that is described as the prototypical 



post-common envelope detached system with a low probability of becoming a Cataclysmic Variable (CV) 



within a Hubble time ( King et al. 1994 and Marks, 1994). This object was selected for our HS T parallax 
program because a directly measured distance could reduce the uncertainty of the radius of one of the 
hottest white dwarfs. Since the instigation of this program and the selection of targets over 15 years ago, 



at least two other groups have measured a parallax for Feige 24; USNO-Flagstaff (Palm et al., 1988|), and 



HIPPARCOS (Perryman et al., 1997 and Vauclair et al., 1997). We outlined the results of a preliminary 



analysis in Benedict et al., 2000. Here we discuss our analysis and final results in detail. 



Provencal et al. 1998 presented radii derived from HIPPARCOS parallaxes for 21 white dwarfs. In 
most cases, the dominating error term for the white dwarf radii was the parallax uncertainty. Our parallax 
of Feige 24, while slow in coming, has provided a fractional parallax uncertainty, — , similar to those in the 



Provencal et al. 1998 study, but for a much hotter, more distant object. 



We time-tag our data with a modified Julian Date, MJD = JD — 2400000.5. We abbreviate millisecond 
of arc, mas; white dwarf, DA; and M dwarf, dM, throughout. 



2. The Astrometry 

Our astrometric observations were obtained with Fine Guidance Sensor 3 (FGS 3), a two-axis, 



white-light interferometer aboard HST. Bradley et al., 1991 provide an overview of the FGS 3 instrument 



and Benedict et al., (1999) describe the astrometric capabilities of FGS 3 and typical data acquisition and 



reduction strategies. 

We use the term 'pickle' to describe the field of regard of the FGS. The instantaneous field of view of 
FGS 3 is a 5 x 5 arcsec square aperture. Figure [l] shows a finding chart for Feige 24 and our astrometric 
reference stars in the FGS 3 pickle as observed on 08 Aug 1997. Note the less than ideal placement of the 
primary science target with respect to the reference frame. The placement of Feige 24, at one side of the 
distribution of reference stars, seems to have produced no adverse astrometric or photometric effects. 
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2.1. The Astrometric Reference Frame 

Table H provides a list of the observation epochs. Our data reduction and calibration procedures 



are described in Benedict et ah, (1999) and McArthur et al., 1999. We obtained a total of 71 successful 



measurements of our reference stars during eight 'observing runs'. For each of these eight observation sets 



we determine the scale and rotation relative to the sky, using a GaussFit ( Jefferys et al., 1988 ) model. The 
orientation of the observation sets is obtained from ground-based astrometry (USNO-A2.0, Monet 1998 ) 
with uncertainties in the field orientation ±0?12. 



Having only 8 observation sets and four reference stars precludes us from our usual practice (|Benedict 



et al., (1999) ) of constraining the proper motions and parallaxes to sum to zero (Y>fi = and S7r = 0) for 
the entire reference frame. From a series of solutions we determined that only reference star ref-3 has a 
statistically significant proper motion and parallax. So, we constrain /i = and n = for reference stars 
ref-2, -4 and -5. 

We conclude from histograms (Figure ||) of the reference star residuals that we have obtained a 
per-observation precision of ~ 1 mas. The resulting reference frame 'catalog' (Table ^) was determined with 
final errors < >= 0.5 and < a, q >= 0.6 mas. 

To determine if there might be unmodeled, but eventually correctable, systematic effects at the 1 
mas level, we plotted the Feige 24 reference frame X and Y residuals against a number of spacecraft, 
instrumental, and astronomical parameters. These included X, Y position within the pickle; radial distance 
from the pickle center; reference star V magnitude and B-V color; and epoch of observation. We saw no 
trends, other than the expected increase in positional uncertainty with reference star magnitude. 



2.2. Modeling the Parallax and Proper Motion of Feige 24 

Spectroscopy of the reference frame stars obtained from the WIYN [] and an estimate of color excess, 
E(B-V), from Burstcin fc Hcilcs 1982| . (Table §) shows that the colors of the reference stars and our science 



target differ, with A(B~ V) ~ —1. Therefore, we apply the differential correction for lateral color discussed 



in Benedict et al., (1999) to the Feige 24 observations and obtain a parallax relative to our reference frame, 



1 The WIYN Observatory is a joint facility of the University of Wisconsin-Madison, Indiana University, 
Yale University, and the National Optical Astronomy Observatories. 
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ir re i — 13.8 ± 0.4 mas. The proper motion relative to the four astrometric reference stars is listed in Table 



Franz et al (1998) and Benedict et al., (1999)| have demonstrated 1 mas astrometric precision for FGS 
3. Table |^ presents our Feige 24 astrometric residuals obtained from the parallax and proper motion model. 
Histograms of these residuals are characterized by a x = 1.0 and a y = 1.2 mas. This was slightly larger 
than expected. To investigate whether or not the larger residuals could be attributed to Feige 24, Figure || 
presents the residuals phased to the VT94 orbital period, P = 4.23160 d , with T = HJD 2448578.3973. We 
find no significant trends in the astrometric residuals. In particular, there is no correlation with the two 



distinct HST orientations required by the pointing constraints discussed in Benedict et al., (1999) . With 
any reasonable masses for the DA and dM components, a binary system at this distance, having this period, 
could exhibit maximum reflex motion at the 0.5 mas level. This null detection does not place very useful 
upper limits on the component masses. 

Because our parallax for Feige 24 is determined with respect to the reference frame stars which 
have their own parallaxes, we must apply a correction from relative to absolute parallax. The WIYN 
spectroscopy and the estimated color excess (See Table 2) indicate a reference frame with an average 
parallax of < tt > r e/ = 0.9 ± 0.4 mas. where the error is based on the dispersion of the individual 
spectrophotometric parallaxes. To check our correction to absolute, we compare it to those used in the 
Yale Parallax Catalog (YPC95, van Altena, Hoffleit, & Lee, Section 3.2). From YPC95, Fig. 2, the Feige 
24 galactic latitude, b = — 50?3 and average magnitude for the reference frame, < V re f >= 13.4, we obtain 
a correction to absolute of 1.9 mas. Rather than use a galactic model-dependent correction, we adopt the 
spectroscopically derived < n > re /= 0.9 ± 0.4 mas. Applying this correction results in an absolute parallax 
of TT a bs = +14.7 ± 0.6 mas, where the error has equal contributions from the HST FGS observations and 
the correction to absolute parallax. Finally, we note that our proper motion is smaller than either the 
HIPPARCOS or USNO values, where the HIPPARCOS value is an absolute proper motion while the USNO 
and the HST values are relative to their respective reference frame proper motions. If our reference stars 
are a representative statistical sample of the parent population, then based on the data in Table III in van 
Altena (1974), we expect a statistical uncertainty in the mean value of the correction to absolute proper 
motion (not applied here) of ±6 mas y -1 . 

We compare our absolute parallax to previous work in Table |] and in Figure [|. We adopt for the 
remainder of this paper the weighted average absolute parallax, < n a b s >= 14.6 ± 0.4 mas, shown as a 
horizontal dashed line in Figure ||. Weights used are I /a 2 . 
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Lutz & Kelker (1973) show that for a uniform distribution of stars, the measured trigonometric 
parallaxes are strongly biased towards the observer (i.e., too large), rendering inferred distances and 
luminosities too small. This bias is proportional to (a^/ir) 2 . Using a space density determined for the CV 
RW Tri (McArthur et al. 1999), and presuming that Feige 24 is a member of that same class of object 
(binaries containing white dwarfs), we determine an LK correction of —0.01 ± 0.01 magnitudes. Correcting 
our distance modulus, we obtain m — M = 4.17 ± 0.11. 



2.3. Kinematic Age of the Feige 24 System 

From the VT94 systemic radial velocity and either our proper motions or those from HIPPARCOS 
(Table |4|) we derive the space velocity of Feige 24, 67 ± 1 km s _1 . The velocity component perpendicular 
to the galactic plane, W, is -37 km s _1 . Our new parallax places the star 53 parsecs below the Sun or 
61 parsecs below the galactic plane. An object this far below the galactic plane and continuing to move 
further away from the plane so swiftly is more characteristic of a 'thick disk' than a thin disk object (c.f. 



Thcjll ct al., 1997 ). Feige 24, if truly a Pop I object, has a space velocity 3.5 times the young disk velocity 
dispersion. These data suggest that Feige 24 formed prior to the formation of the galactic disk, although 
subsequent evolution of the DA component is likely quite recent. This may be an instance of past mass 
transfer in an intermediate Population II object. 



3. Astrophysics of the Feige 24 System 

We discuss the consequences of a more precisely determined parallax, calculating some astrophysically 
relevant parameters for the DA and dM components. These are collected in Table |[ Our goals are the 
radius and mass of the DA component. We first calculate a radius, then estimate the time since the DA 
formation event. Component masses have been estimated by VT94. We will revisit this issue later. That we 
do not substantially improve the mass uncertainty motivates a future direct measurement of the component 
separation. This one measurement would yield precise masses. A series of measurements would provide 
individual orbits, possibly illuminating past and future component interactions. 
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3.1. Estimating the DA Radius 



To estimate the DA radius we require an intrinsic luminosity. From Landolt 1983 we obtain a system 



total magnitude, Vtot — 12.41 ± 0.01. The magnitude of the white dwarf is critical and difficult to obtain, 



because the M dwarf always contributes flux. Holberg et al., 1986 derive Vda — 12.56 ± 0.05 using IUE 



spectra. They ratio Feige 24 with other hot DA; G191 B2B, GD246 and HZ43. From the DA magnitude 
and total magnitude we obtain V^m = 14.63 ± 0.05 and AV^ = 2.07. We assume an Ay = for Feige 24 at d 
= 69 pc, consistent with our adopted Ay = 0.09 for the reference frame at an average distance d = 1600 pc 
(Table |2|). The LK bias-corrected distance modulus [m — M = 4.17 ±0.11) then yields absolute magnitudes 
My = 10.46 ± 0.12 for the red dwarf companion and My = 8.39 ± 0.12 for the DA. 

A recently determined temperature of the Feige 24 DA, taking into account non-LTE and heavy 
element effects ( Barstow Hubeny fc Holberg 1998 ), is T®A = 56, 370 ± 1, 000K. This temperature yields a 



radius via differential comparison with the sun. This procedure requires a bolometric magnitude, hence, 
a bolometric correction. We could adopt the bolometric correction, B.C. = -4.88, generated by Bergeron 
et al. (1995) from a pure Hydrogen, log g = 8 DA model convolved with a V bandpass. But, Feige 24 is 
neither log g = 8 nor pure H. 

Flower (1996) provides bolometric corrections for normal stars up to T e ff ~ 54,000^. From Flower 
(1996), figure 4, the relationship between logT e ff and B.C. is linear for T e ff > 25,000-ftT. Hotter stars lie 
on the Rayleigh- Jeans tail of the blackbody curve, where flux is roughly proportional to T e ff, not T^fj- A 
small linear extrapolation yields B.C. = -4.82±0.06 for the Feige 24 DA. The B.C. error comes from the 
uncertainty in Because a DA with some heavy elements in its atmosphere radiates more like a hot 

normal star than a pure H DA, we choose the Flower correction rather than the model correction. We are 
also encouraged by the near equality of the B.C. from observation and theory. 

We obtain a DA bolometric luminosity M^f = My + B.C. = 3.57 ± 0.13. R DA follows from the 
expression 

M? ol - M h D Q f = 10 log(T%/ T§ f ) + 5 log(R DA /R Q ) (1) 
where we assume for the Sun A/® = +4.75 and T® ff = 5800K. We find R DA = 0.0180 ± O.OO13i? , 



following the error analysis of Provencal et al. 1998 . The primary sources of error for this radius are the 



bolometric correction and the T®A . 

A second approach to deriving Rda involves the V-band average flux, Hy, discussed in Bergeron et al. 
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(1995). They list H$ A as a function of temperature for, again, the pure Hydrogen, log g = 8 model. If we 
can determine an H ® , we can derive Rda from 

R DA = (H^/H^lO-^r-M®) (2) 

where M$ A = 8.39±0.12 comes from our parallax and My = 4.82 is assumed. We obtain Hy by convolving 
the Bessel (1990) V band response with the solar spectral distribution listed in Allen (1973). We calculate 
if© = 6.771 x 10 5 ergs cm" 2 s" 1 A" 1 str" 1 . We obtain for T eff = 56,370^ an R DA = 0.0188 ± 0.0010. 

A weighted average of the two independent determinations provides Rda = 0.0185 ± 0.0008i?Q, where 
the error is certainly underestimated due to unknown systematic effects. Parallax is no longer a significant 
source of error for the radius determination. Comparing with the results presented in Provencal et al. 1998j , 
figure 7, we find Feige 24 to have a radius larger than any other white dwarf. 

With a temperature T e // ~ 56, 000, the time since the DA formation event is unlikely to be longer than 
1.5 My. This conclusion is drawn from the DA cooling tracks as function of mass calculated by M. Wood, 
detailed in Sion 1999 , fig. 7. These models also indicate that the DA mass must satisfy A4 da > 0.47Wq to 
remain near this lofty T e f f for longer than 3 x 10 5 y. 



3.2. Estimating The White Dwarf Mass 

Before estimating A4 da we review the VT94 minimum component masses from their radial velocities 
and the Kepler relation for total system mass, separation, and period. Then, we estimate the DA mass 
using two different approaches. We first attempt to determine the most likely dM mass. The VT94 radial 
velocity amplitude ratio then provides the DA mass. The second, independent mass estimate follows from 
our derived radius along with the DA atmospheric parameter, log g, obtained through spectroscopy. Our 
DA mass estimate will differ little from VT94, and, if better, is so only by virtue of more recent dM models 
and DA atmospheric parameters. 



3.2.1. Minimum Component Masses from Binary Radial Velocities 

The system total lower mass limit can be set by the VT94 radial velocities and the Kepler relation 
for mass, separation, and period. VT 94 give us the velocities along each component orbit, the fact 
that each orbit is circular (from the pure sine wave fits to the velocity curves), and the period, the time 
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it takes to travel around each orbit. Assuming an edge-on system (i = 90deg), one that can produce 
the full vector amount of radial velocity amplitude measured by VT94, the minimum system mass is 
Mtot = 0.73.Mq. From the VT94 mass ratio, McLm/Mda = 0.63 ± 0.04, we obtain the DA mass limit, 
■M-da Kepler > 0.44Mq, and the dM mass limit, MdM Kepler > 0.26.Mq. No smaller masses can produce 
the observed radial velocities for orbits of these known sizes. At d = 68.5 pc an edge-on system with 
minimum mass would separate the components by 672 microarcsec, or 9.9Rq. 



3.2.2. Inclination from the Light Curve 

VT 94 find Ha equivalent width variations that phase with the orbital period. These show a maximum 
at cj> = 0.5. Photometric variations of Feige 24 might be detectable, because the photometric capabilities 
of FGS 3 approach a precision of 0.002 magnitude (Benedict et al., 1998). Figure [| shows the flat-fielded 



counts and the corresponding differential instrumental magnitudes as well as a sin wave fit with amplitude 
and phase as free parameters. There is a clear photometric signature with a peak-to-peak amplitude 0.028 
magnitude, showing maximum system brightness at phase <j) — 0.58 ± 0.09. Given the sparse coverage, this 
phase at maximum is not surprisingly different from the Ha equivalent width maximum seen at (j) = 0.5. 

A likely mechanism for producing the single-peaked orbital light curve is heating of the dM star by 
the white dwarf (the reflection effect). As the dM star orbits the white dwarf, its heated face is alternately 
more or less visible, increasing and decreasing the observed flux from Feige 24 once per orbit. To test this 
hypothesis we calculated model light curves using an updated version of the light curve synthesis program 
described by Zhang et al. (1986). We initially adopted T eff = 56, 370 K and R — 0.0185 R Q for the white 
dwarf, T eff = 3800 K and R — 0.52 R Q for the M1-2V star, and 4.8 x 10~ 2 AU (10.3i? Q ) for the separation 
of their centers of mass, and then adjusted the temperature of the dM star so that it contributed 13.5% of 
the V flux from the system. The peak-to-peak amplitudes of the resulting model light curves are a function 
of orbital inclination, topping out at ~ 0.025 mag for i = 90°, and can easily be made to agree in amplitude 
and shape with the observed light curve. 

This photometric behavior is entirely consistent with reflection effects (c.f. Robinson et al. 2000). 
We find that the quality of the observed light curve is, however, inadequate to improve the parameters of 
the system, particularly the inclination. We have not sufficiently sampled the expected flat section of the 
light curve (near = 0). Nevertheless, these results do provide quantitative evidence that (1) the orbital 
light curve is caused by heating and (2) the heating is consistent with the radius and temperature we have 
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derived for the white dwarf - a useful external check on our results. 



3.2.3. DA Mass from the M Dwarf 



The dM absolute magnitude (My = 10.46 ±0.12 ) implies a spectral type M2V (Henry, Kirkpatrick, & 



Simons 1994), consistent with Liebert & Margon 1977. The absolute magnitude of an M dwarf star depends 



not only on mass, but also on age (evolutionary stage) and chemical composition. Baraffe et al., 1998 have 
produced a grid of models, varying metallicity, [M/H] and helium abundance, Y. We plot in Figure ^ their 
Mass-Luminosity curves for dwarfs with ages lOMy and lOGy, with [M/H] = and lOGy with [M/H] = 
-0.5, all with solar helium abundance. The complete grid of Baraffe et al. models shows that M dwarfs in 
the mass range 0.175 < MdM < 0.43A4q with -0.5 < [M/H] < have My = 10.46 at some time in their 
evolution from lOMy to lOGy. 

The dM mass now depends on metallicity and how quickly an M dwarf of a given mass decreases 
in brightness. Figure shows the dependence of brightness on mass, age, and metallicity. These Baraffe 
et al. models indicate that solar metallicity stars with higher mass remain near My = 10.46 far longer 
than low mass stars. However, kinematically, Feige 24 is more likely to be old and of lower than solar 
metallicity than young and of normal metallicity. First adopting the lOGy model, [M/H] = 0, and calculated 
absolute magnitude, we estimate the dM star mass MdM = 0.43 ± 0.08.Mq, because that mass remains at 
My = 10.46 for a larger fraction of the total lifetime than any other. However, if we accept the kinematical 
suggestion of allegiance to a thick disk population, then [M/H] < is more likely. Assuming [M/H] = -0.5 
results in a dM star mass MdM — 0.185 ± 0.08A4q. 

Radial velocities from VT94 (dM from Kitt Peak, DA from IUE) provide the velocity amplitude ratio, 
K DA /K d M = 0.63±0.04 = MdM / Mda- From the total possible dM mass range, 0.185 < MdM < 0.43.Mq, 
and the mass ratio we derive a DA mass range, 0.29 < Mda < 0.68Mq. Applying the limit, 
Mda Kepler > OAAMq, we obtain 0.44 < Mda < 0.68A4Q. Keplerian lower limits argue for a dM star 
mass 0.26 < MdM < 0.43.Mq, a range consistent with a metallicity slightly less than solar and an age in 
excess of 0.3 Gy (Figure 0). 
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3.2.4- DA Mass from Atmospheric Parameters 

The dM star does not provide a particularly precise DA mass estimate. If one knows the surface 
gravity, g, and the radius, R, the mass can be obtained through 



M = gR 2 /G 



(3) 



where G is the gravitational constant. The quantity log g comes from analysis of the line profiles in spectra. 



Recent determinations include: Marsh et al. 1997, log g = 7.53 ± 0.09; Kidder 1991, log g = 7.45 ± 0.51 



Vcnnes et al., 1997, log g = 7.2 ± 0.07; Finley Koester fc Basri, 1997, log g = 7.17 ± 0.15; and Barstow 



Hubcny fc Holberg 1998) , log g = 7.36 ± 0.12. The full range of the measures and Equation 3 yield the range 



of mass values 0.21 < Mda ^ 0.47A4q. Applying the limit, AAda Kepler > 0.44A4q eliminates nearly all 
of these mass determinations. In this case our radius and the Kepler limit indicate that log g should be at 
the high end of these measures. 



3.3. The White Dwarf Composition 

We next place Feige 24 on a white dwarf mass-radius diagram (Figure ||) . We plot our two independently 
determined mass ranges against our adopted radius, Rda — 0.0185 ± 0.0008i?Q. We represent the radius 
error by the two horizontal long-short dashed lines. The top thick horizontal bar shows the M. da 
determined from atmospheric parameters. Only the largest log g at the largest radius produces masses 
in excess of the Keplerian limit. The thick bar at Rda — 0.0185i?Q indicates the Mda range derived 
through the dM mass estimates. For this determination the mass error bars indicate the range of ages and 
[M/H] discussed in section 3.2.2. For any dM older than 1-2 Gy the lower masses are associated with lower 
metallicity. The vertical bold dotted line shows the lowest possible Mda that can produce the observed 
VT 94 radial velocity amplitudes for an edge-on orientation of this binary system. We also plot several 
values of log g (dashed) and j gra v (thin solid) . The wide grey curves in Figure || are C and He DA models 



from Venncs et al., 1995. While uncertain, a carbon core DA seems more likely than a pure He core DA. 



4. Discussion 



While our estimated dM and DA masses differ little from VT94, our DA radius differs substantially. 
VT94 note the difference between their minimum radius, Rda — 0.028-Rq, and that predicted by the Dahn 
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et al., 1988 parallax. This discrepancy is exacerbated by the two new parallax determinations (HST and 
HIPPARCOS), folded into our weighted average parallax. 

VT94 derive a DA gravitational redshift, j gra v = 8.7 ± 2 km s _1 from the measured mean velocities 
for the dM and DA. Combined with our Rda = 0.0185i?Q, this 7 sro u suggests a forbidden DA mass, 
Mda ~ O.3A^0. Reducing the mass of the DA component could reconcile the VT94 log g and jgrav with 
our radius. 

We speculate that a third component in the Feige 24 system, a low-mass companion to the DA star, 
could preserve the total system mass and lower the DA mass. If all components are coplanar, the VT94 
DA radial velocities apply strict limits to this reconciliation, because too high a mass for component C 
would show up as large residuals. We estimate from the scatter that a radial velocity amplitude of ±10 
km s _1 could 'hide' in the VT94 DA radial velocity measurements. Stellar dynamics applies yet another 
constraint. Holman & Wiegert (1999) parameterize the stability of tertiary companions as a function of 
stellar component A and B mass function, /j, = Ma/ (Ma + -Mb), and AB binary orbit ellipticity, e. With 
c=0 and /i = 0.39 we find (from their table 3) that component C must have an orbital semi-major axis less 
than 0.3 times that of AB. 

Insisting that Mda — O.3O.M0 (this mass - with our radius - would produce the upper limit VT94 
Igrav = 10.7 km s _1 ) requires Mc = 0.14M© {Ma + Mc = O.MMq). To hide the C component from 
the radial velocity technique requires a very low AC inclination, nearly face-on. However, non-coplanarity 



reduces the size of the stable AC semi-major axis even further (Weigert & Holman, 1997; Pendleton & Black 



1983| ). As an example suppose component C must have an orbital semi-major axis of 0.1 or less that of AB 



to insure stability. An AC period, P = 0.18 d (A.3 h ), and i = 6° would produce a radial velocity signature of 
about ±10 km s _1 . Finally, the Mass-Luminosity Relation of Henry et al. (1999) would predict My = 14.0, 
hence, Vc ~ 17.2, likely undetectable in any of the spectra analyzed for radial velocities. Have we built a 
new CV, one that should evidence mass transfer and all the associated phenomena? A recent review of CVs 



( |Bcucrmann, 1999| ) indicates that the putative component C (Mc = 0A4Mq) would have to orbit much 
closer (P~ l.h h ) to the DA primary before filling its Roche lobe and producing the characteristic signature 
of a CV. 

Finally we note that our radius differs little from that derived by VT94 from the only trigonometric 



parallax then available (Dahn et al., 1988). The unresolved inconsistency between radii (derived from direct 



parallaxes) and surface gravities (derived from minimum mass and those radii) illuminates the need for 
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high angular resolution observations and direct mass determinations. 

The Feige 24 DA mass will rest on an age- and metallicity-dependent lower main sequence Mass- 
Luminosity Relationship or still uncertain log g measurements until the component separations are 
measured directly. Resolving the inconsistencies between the DA mass estimates (involving dM stellar 
models and uncertain temperatures, log g, and bolometric corrections) requires astrometry, both to further 
reduce the parallax uncertainty, and, more importantly, to spatially resolve this system. Astromctrically 
derived orbital parameters will provide unambiguous and precise mass determinations for both components. 
They may also offer insight regarding past and future component interactions. 

This system and dozens more like it are ideal targets for the Space Interferometry Mission 



( http: / / sim.jpl.nasa.gov ) . Feige 24, at a distance of 69 pc with P = 4.23 d , has a total component separation 
on order 700 microarcsec. The component orbits are much larger than the expected SIM measurement 
limits. Because shortward of 700 nm 70-80% of the system flux is contributed by the DA (Thorstensen et al. 
1978), the wide SIM bandpass and spectral resolution should allow measurement of positions, magnitudes, 
and colors for both components, even with AV^ ~ 2. 

Once launched SIM will provide crucial astrometry for this and similar systems at ten times the 
distance (determined by target magnitude, not astrometric precision). SIM measurements of this system 
along with many other binaries will provide data with which to create an age- and metallicity-dependent 
Mass-Luminosity Relationship of exquisite accuracy. 



5. Conclusions 

1 . The weighted average of three independent parallax mesurements yields a distance to the dM + DA 
binary Feige 24 with n D /D = 2.8%. D = 68.4±^;° pc. 

2. We estimate the radius of the DA component using two methods. The first requires either a 
model-dependent bolometric correction, or one that derives from hot, normal stars. The second utilizes a 
model-dependent V-band average flux, Hy. The two results agree within their errors and yield a weighted 
average Rda = 0.0185 ± 0.0008i?Q, where the most significant contributions to the error are the uncertain 

L eff 



TR4 and B.C.. This radius is larger than any of the WD discussed in Provencal et al. 199S 



3. FGS photometry provides quantitative evidence that the orbital light curve is caused by heating of 
the dM component by the DA. That signature is consistent with the assumed temperature and the radius 
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we have derived for the white dwarf 

4. The VT94 measured radial velocity amplitudes, amplitude ratios, and the assumption of Keplerian 
circular motion exclude M.da Kepler < 0.44.Mq and Kepler < 0.26.Mq. 



5. We estimate the dM component mass, 0.26 < AidM < 0.43.MQ, from the Baraffe et al., 199S stellar 
evolution models, a lower limit from Keplerian circular orbits, and the VT 94 radial velocities. The upper 
range is due to unknown age and metallicity, [M/H]. A DA mass range (0.44 < M.da < 0.68.Mq) follows 
directly from the VT94 radial velocity amplitudes. 

6. We determine A4da from our Rda and a rather wide range of spectroscopically determined 
log g values. This approach yields 0.44 < Mda < 0.47A4q, where again the lower limit is imposed by 

MdA Kepler > OAAMq. 

7. We plot these DA component mass ranges on the M. — R plane. With the assistance of the hard 



lower mass limit and C and He DA M. — R models from Vennes et al., 1995 , we identify Feige 24 to have a 
carbon core. A pure He core DA seems less likely. 

8. Noting that our radius and the minimum possible M.da are inconsistent with the VT94 j gr av, we 
explore the possibility of a tertiary component. A component C, orbiting a common center of mass with 
the DA, having a period in the range 1.5 <P < 5 h with the orbit plane nearly face-on, could reduce the DA 
mass to M.da = O.SOA^q and not produce any observational evidence. 

8. SIM will be able to measure the orbits of each known component and provide directly measured 
dynamical masses for both. Orbit size and precise shape may provide information on the nature of past and 
future interactions between the two components. SIM would also detect a tertiary, if present. 
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Table 1. Feige 24 Reference Frame: Astrometry 



ID 


V 




n 




/iy 


7T 






(arcsec) 


(arcsec) 


(arcsec y^ 1 ) 


(arcsec y _1 ) 


(arcsec) 


ref-2* 


11.59 


0.0±0.0004 


0.0±0.0004 











ref-3 


13.38 


-118.9943±0.0004 


50.4968±0.0004 


0.0157±0.0004 


- 0.0014±0.00050 


-0.0006±0.0003 


ref-4 


14.82 


27.3168±0.0010 


63.5709±0.0010 


- 








ref-5 


13.66 


-144.4703±0.0004 


86.9851±0.0005 












*RA, Dec = 288.087967, 2.898281 (J2000) 



Table 2. Feige 24 and its Reference Frame: Stellar Parameters 



ID 


V a 


B-V h 


SpT 


My 


E(B-V) c 


Ay 


m-M 


D(pc) 


n abs (mas) 


ref-2 


11.59 


1.00 


G9 III 


0.75 


0.03 


0.093 


10.84 


960 


0.6 


ref-3 


13.38 


0.69 


G3 V 


4.8 


0.03 


0.093 


8.58 


410 


1.9 


ref-4 


14.82 


0.61 


GO V 


4.4 


0.03 


0.093 


10.42 


935 


0.8 


ref-5 


13.66 


0.63 


F9 III 


1.2 


0.03 


0.093 


12.46 


2540 


0.3 


Feige 24 


12.41 d 


-0.20 d 

















a from FGS PMT measures calibrated as per Nelan et al.1999 
b from B-V= f(Sp.T.) + E(B-V) 



c from Burstcin & Hciles 1982 



from Landolt 1983] 
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Table 3. HST Observations of Feige 24 and Astrometric Residuals 



Obs. Set MJD X residual Y residual 



1 


49930.92188 


0.0002 


0.0008 


1 


49930.9375 


0.0000 


-0.0012 


1 


49930.94531 


-0.0003 


-0.0004 


2 


49936.88672 


0.0005 


-0.0002 


2 


49936.90234 


-0.0002 


-0.0003 


2 


49936.91016 


0.0006 


0.0001 


3 


50102.09375 


-0.0014 


0.0008 


3 


50102.10938 


0.0001 


-0.0002 


3 


50102.11719 


0.0004 


-0.0013 


4 


50109.06641 


-0.0008 


0.0007 


4 


50109.07813 


-0.0001 


0.0002 


4 


50109.08594 


0.0007 


0.0001 


5 


50669.78125 


-0.0001 


0.0010 


5 


50669.79688 


0.0009 


0.0016 


5 


50669.80469 


-0.0010 


-0.0004 


6 


50678.92188 


-0.0003 


-0.0004 


6 


50678.9375 


0.0003 


0.0006 


6 


50678.94531 


-0.0016 


0.0010 


7 


50819.28125 


0.0014 


-0.0003 


7 


50819.28906 


0.0005 


-0.0008 


7 


50819.29688 


0.0003 


0.0001 


8 


50821.29688 


0.0006 


-0.0009 


8 


50821.30469 


-0.0005 


-0.0008 


8 


50821.3125 


-0.0007 


0.0006 
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Table 4. Feige 24 Parallax, Proper Motion, and Radial Velocity 



Parameter 



Value 



HST study duration 

number of observation sets 

Number of ref. stars 

ref. stars < V > 

ref. stars < B - V > 

HST Relative Parallax 

corr to absolute 

HST Absolute Parallax 

HIPPARCOS Absolute Parallax 

USNO Absolute Parallax 

HST Proper Motion (p) 

in p. a. 

HIPPARCOS ft 
in p. a. 
USNO n 
in p. a. 



2.4 y 

8 
4 

13.4 
0.7 

13.8 ± 0.4 mas 
0.9 ± 0.7 mas 

14.7 ± 0.6 mas 

13.4 ±3.6 mas 

13.5 ± 2.9 mas 
71.1 ± 0.6 mas y" 1 
83?6 

85.8 ± 5 mas y" 1 
84? 2 

78.4 ± 1.9 mas y" 1 
88?4 



Weighted Average Absolute Parallax 14.6 ± 0.4 mas 

m-M (LK bias corrected) 4.17 ± 0.11 

System Radial Velocity, 7 + 62.0 ± 1.4 km s" 1 VT94 

Galactocentric z velocity, W -37 ± 1.5 km s" 1 7 & HST or HIPPARCOS n 
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Table 5. Feige 24 Astrophysical Quantities 



Parameter 



Value 



Source 



Vtot 

B-V 
V DA 

VdM 

A v 

m-M (LK bias corrected) 
dM M v 
dM Sp. T. 

M dM 
DA M v 
DA B.C. 

M DA 

m bol 

rpDA 

1 eff 

Rda 

MdM IM D a 

Mda 

Mda 



12.41±0.01 
-0.20 ± 0.01 
12.56 ±0.05 
14.63 ±0.05 


4.17±0.11 
10.46 ±0.12 
M2V 

0.29 - O.437U 
8.39 ±0.12 
-4.82 ±0.06 
3.57 ±0.13 
56,370±l,000if 
0.0185 ± O.OOO8i? 
0.63 ±0.04 
0A9t° o f 5 M @ 
0.44 - O.477W 



Landolt 1983 



Landolt 1983 



Holberg et al., 1986| 



Vtot & Vda 

reference frame < B — V > 

this paper 

m-M 



Sp.T. (Table |) 



dM M v & Henry, Kirkpatrick, fc Simons 1994 



dM M v , [Baraffc ct al., 1998| , M da Kepie 
m-M 

Flower, 1996 

= DA M v + B.C. 



Barstow Hubeny & Holberg 1998 



this paper 

K DA /K dM VT94 

MdM, K D aI 'KdM& Mda Ke P u 

log g & Mda Kepler 
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File : OzTeX: to-from cfgauss, abridged:Astrophys ji:Feige 24 7i:feige coFile : Yzeg>A Pickles 
V1: Ra: 139731.543" Dec: 12659.112" Roll: 13.88° Orient: -92.71° Veh.Roll: -87.29° 
AntiSun: 80° Moon: 164° Plate Roll: -1° Tobs: 97/08/09.500 Teat: 00/01/01.500 Now: 99/7/12 13:41:01 




Fig. 1. — Location of reference stars within the FGS 3 field of regard on 08 Aug 1997. Note the less than 
ideal placement of the primary science target with respect to the reference frame. 
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Y Residuals 
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Fig. 2. — Histograms of x and y residuals obtained from modeling the Feige 24 reference frame to obtain 
scale, orientation, and offset parameters. Distributions are fit with gaussians. 
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orbital phase (P=4.23160 ) 



Fig. 3. — Astrometric residuals in RA (X) and Dec (Y) phased to the VT94 orbit (P = 4.23160 d and 
To = JD2448578.3973). Boxes and open circles denote the two HST orientations, which seem to have no 
effect on the astrometric residuals. Dashed lines are best-fit sine waves constrained to the VT94 period. 
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<7i abs > = 14.6 ± 0.4 mas 



Fig. 4.— Absolute parallax determinations for Feige 24. We compare HST, HIPPARCOS, and USNO (Palm 



ct al., 1988| ). Error bars are 1-a. The horizontal dashed line gives the weighted average absolute parallax, 

< Hubs >■ 
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orbital phase (P=4.23160 d ) 



1.2 



Fig. 5. — Flat fielded intensity (the filter (F583W) has a bandpass centered on 583 nm, with 234 
nm FWHM) and differential instrumental magnitudes phased to the VT94 orbit (P = 4.23160 d and 
To = JD2448578.3973). Boxes and open circles denote the two HST orientations, which seem to have 
no effect on the photometry. The dashed line is a best-fit sine wave constrained to the VT94 period. 
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Fig. 6. — M dwarf absolute magnitude as a function of mass from the stellar evolution models of Baraffe 



ct al., 1998| . A wide range of masses, ages, and metallicities can result in the derived dM absolute magnitude, 
Mv — 10.46 (horizontal dashed line). Note that at a given mass, a low metallicity ([M/H]=-0.5) star is always 
brighter than a high metallicity ([M/H] = 0.0) star. 
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Fig. 7. — The time variation of absolute magnitude for M dwarfs of various masses taken from the stellar 



evolution models of Baraffe et al., 1998. Empty symbols are for solar metallicity; filled for low metallicity 
([M/H]=-0.5). Note that the higher mass stars remain near the computed M dwarf absolute magnitude, 
My = 10.46 (dashed line), far longer than the low mass stars. A low metallicity star with AidM — 0.175.Mq 
is brighter than a solar metallicity star with AidM = OAOMq. 
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Fig. 8. — Feige 24 mass and radius on a DA mass-radius map. The vertical bold dotted line shows the lowest 
possible M. da from keplerian considerations. The radius error is represented by the top and bottom long- 
short dashed lines. The top thick horizontal bar shows the M.da determined from atmospheric parameters. 
Only at the largest radius (lowest temperature) and largest log g do we obtain a DA mass in excess of the 
keplerian limit. The bottom thick horizontal bar at Rda = 0.0185i?Q indicates the Mda range derived 
through the dM. We also plot several values of log g (dashed) and j grav (thin solid). The grey wide solid 
lines are C and He DA models from IVcnncs ct al., 1995[ A C core DA is somewhat more likely than an He 
core DA. 



